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Abstract 

The effect of hard processes to be encountered in HBT studies at the 
Large Hadron Collider have been studied. A simple simulation has allowed 
us to generate momentum correlations involving jet particles as well as 
particles originating from the kinetic freeze out and to compare them to 
a simple theoretical model which has been developed. The first results on 
the effect of hard processes on the correlation function for the case of jet 
quenching are presented. 



1 Introduction 

Momentum correlations of pions and kaons created in heavy ion collisions have 
been extensively studied in a wide energy range from the AGS through SPS and 
finally at RHIC pQ. The purpose of these studies is to get information on the 
dimensions of the phase space from which the particles have been emitted |2J. 
However, interpretation of the data measured at RHIC has been proven to be 
problematic and is not yet completely understood (see e.g. 
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Looking ahead towards LHC one can wonder what are momentum correlations 
going to look like in presence of an important contribution of hard processes. If 
we assume that the source is a superposition of jet sources and a thermal fireball 
then the resulting correlations should reflect contributions from: 

• pairs of particles from a single jet, which will reflect dimension of the region 
where the jet fragmented; 

• pairs of particles from different jets, which will be given by the size of the 
initial collision volume; 

• pairs where both particles stemming from the decoupling of thermal fireball 
(larger than the initial collision volume); 

• pairs where one particle comes from the thermal fireball and the other one 
is produced from jet fragmentation. 

In this paper we describe a toy model to simulate the effects listed above. We 
provide analytical calculations of the correlation function in a special limiting 
case. Then we investigate how particle production from jet fragmentation shows 
itself in the correlation data. We study the visibility of the effect at different 
proportions of jet production to the total multiplicity and at different transverse 
momenta. 

2 Simulation 

Jets. With Pythia |4 we simulate p+p collisions at y/s = 5.5 TeV. We extract 
jets in pseudorapidity coverage of the ALICE detector — 1 < rj < 1, by making 
use of a simple clustering algorithm. We identify the highest pt particle among 
fragmentation products of a single string as the leading one and then search 
for particles with an opening angle A<p and relative rapidity Arj to the leading 
one that fulfills criterion: a/ A<fi 2 + Arj 2 < 2. Jets with the sum of ptS from all 
particles between 5 and 10 GeV are selected for further study. 

Within the jet, the original position of the particles is distributed accord- 
ing to a Gaussian of the width J and all particles are assumed to be produced 
instantaneously 
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The jet itself is placed into a random position (xj, tj). 
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Jet quenching model. We model the strong jet suppression by a very 
dense medium by assuming that the jets originate from surface of a cylinder with 
radius R and length L. Eventually, we allow for some finite thickness 5R for 
the cylinder surface. It is assumed that the total transverse momentum of the 
jet is always perpendicular to the cylinder surface and that jets are produced 
inst ant aneously. 



Non-quenching model. If the early medium does not suppress jets, we 
should be able to observe jets from hard collisions in all reaction volume. We 
model such a case by choosing a Gaussian distribution of the jets in the transverse 
plane, but we keep the uniform distribution in longitudinal direction from —L/2 
to L/2. Instantaneous production is assumed. 



Background. Apart from jet fragmentation, particles can be produced by ther- 
mal fireball. If we are interested in jets, these particles represent background for 
us. Its momentum distribution is such that the total resulting spectrum at low 
p t is given by 

\ d 2 N ( m t \ 

F(Pi) = ^^ KeXP v-^J' (2a) 

where T is an effective slope parameter. The background spectrum is defined 
as F b (p t ) = F{p t ) — Fj(p t ), where Fj{p t ) the spectrum produced by jets. 

At high p t the spectrum is dominantly shaped by jet fragmentation. Tech- 
nically, it eventually shows up when Fj(p t ) supersedes parametrisation (f2~aj) . In 
order to obtain a smooth total spectrum with exponential low-p t part and high-p t 
tail as given by jets, we choose the normalization of the background spectrum 
Fbipt) = F(p t )/2 when the jet contribution becomes Fj{p t ) > F{p t )/2. 

The background particles are assumed to be produced instantaneously at the 
same time as the jets and are distributed in space according to a Gaussian 



Sfap) = F ^ \ 27lB 2 {pt)) m «p ) 6{t - to) ■ (2b) 
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The size of the background source first decreases linearly with p t up to pt = p\ 
and than stays constant: 
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r B ^^ + for p t <pf^ 



^ ^max* for Pt >p? aX ' 



where we can specify B , -B max , and pf 



Correlation function. These are constructed in the so-called out- side-long 
coordinate frame by making use of the weighting algorithm due to Lednicky j3] 
within the ALICE HBT analysis package HBTAN 0. 
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3 Theoretical understanding 



If we specify the emission function of the source S(x,p), correlation function can 
be calculated from 

I f d A x exp(z q ■ x) S(x, K) | 2 . , 

C(q,K)-l = y- — 1 , n . 3 

If there are two or more jets in an event, and we sum over a large number of events, 
we can express the emission function by integrating over all possible locations of 
the "jet centres" 

S( X , P )«Jd Xi D( Xj , P ) s(x ; Xj ), (4) 

where s(x;xj) was specified in eq. ((H) and D(xj,p) is the distribution of jet 
centres, which lead to production of momentum p. There may be additional 
dependence of S(x,p) on the momentum, but this is irrelevant for calculating the 
correlation function. 

For the case of jets produced at cylinder surface we write 

D( Xj ,p) = l){i,.r r p ,.:.,) = 

exp (^jp) Sfa - R) e( Zj - L/2) Q(L/2 - z 3 ) Sfe - t ) . (5) 

For simplicity, we neglect the possibility of finite thickness of the surface here. 
The factor exp(cos (pj/$) allows particles to be produced under some angle with 
respect to the jet angle. Since we assume that the jet momentum is perpendicular 
to the surface, and we work in the out-side-long system, ipj is the angle between 
Pt of the particle and the radial jet coordinate r. Parameter $ regulates the 
broadness of the jet: $ < 1 leads to very focussed jets, while $ — > oo leaves 
us with uncollimated particle production. In principle, $ depends on particle 
momentum, but we will suppress this dependence here for simplicity. 
Correlation function can be calculated from eq. (JSJ): 

C(q , q s ,qi)-l = exp (-(ql + q 2 s + qf)J 2 ) sin 2 



/7T 
dcp exp [i R (q Q cos <p + q s sin cp)] exp 
-7T 



(2^(1/$) 

2 



(6) 



where Io is the modified Bessel function. For qualitative understanding it is useful 
to examine the limit $ — > oo when one can analytically calculate the cuts of the 
correlation function along the g-axes 

C(q ,q s = 0,qi = 0) - 1 = exp(-J 2 g 2 ) { J 2 (g G R) + H 2 (|g Q R\) } , (7a) 
C(q o = 0,q s , qi = 0)-l = ex V (-J 2 q 2 s )Jl(q s R), (7b) 



C{q Q = 0,q s = 0, q{) - 1 = exp(- J 2 q 2 ) -^j^ sin 2 ( ^ ) , (7c) 
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Figure 1: Correlation function resulting from simulation of 10 jets produced at 
the surface of a cylinder with the radius 6 fm and length 1 fm. Plotted are cuts 
along g-axes. 



where Jo is a Bessel function and H is a Struve function. 

If a background source is added, the correlation function includes contribu- 
tions of pairs of jet particles, background particles and those with one particle 
from a jet and one from background. The analytical expression thus becomes 
more complicated. 

Finally, we also made simulations with jets produced from a Gaussian distri- 
bution. Calculations similar to those in this section show that such a distribution 
of jets leads to a Gaussian correlation function. 



4 Results 

Only jets. First, we analyze events that contain jets only (and no background). 
We assume cylindrical geometry with radius of R = 6 fm, 5R = 1 fm, and length 
L = 1 fm (Fig. [TJ) or L = 2 fm (Fig. |2J). Jet radius J is 1 fm. 

In Figures n and 121 we show the projections of the correlation functions from 
simulations with 10 and 100 jets, respectively. In the latter case the simulated 
correlation function is rather well reproduced by theoretical curves which were 
calculated from eq. (pj). A parameter A < 1 which multiplies the theoretical 
expression was added in order to account for the intercept of C(q — 0) which 
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Figure 2: Correlation function resulting from simulation of 100 jets from a surface 
of a cylinder with the radius 6 fm and length 2 fm. Curves correspond to a fit 
with the theoretical curve given by eq. © multiplied by interception parameter 
A = 0.77. Parameters of the curves are: R = 5.92 fm, J = 1.06 fm, L = 1.82 fm, 
and $ = 0.86. 



is smaller than 2 because in obtaining the projections we always integrated over 
finite regions in the remaining two g-coordinates. 

The structure of the projection in q s in Fig. |21 can be understood from the 
$ — > oo limit formula, eq. (j7b|) : the sub- leading peaks stem from the Bessel 
function. A similar structure in q Q is "dissolved" faster when the jet collimation 
is finite and leads to the "shoulder" in Fig. El (left). Projection in q\ can be 
reproduced by eq. (fTcj) . 

In the simulation, the angular width of the jet depends on transverse mo- 
mentum of particles. Thus $ is not constant. It appears that averaging over 
&(pt) is more important in case of a smaller number of jets. Projections in q Q 
and q s in Fig. Q result from such averaging and we cannot describe them with 
the prescription ©. Beyond the leading peak we do not see any clear subleading 
maxima, but the correlation function assumes a shoulder-like shape at larger g's. 

Jets with background. If in addition to jets there is also other source of 
particles, the characteristic features of the correlation function due to jets get 
diluted. We study this by adding a background source, as described by eqs. (J2J). 
We choose T = 210 MeV, B = 10 fm, 5 max = 1 fm, and p t max = 1.2 GeV. 
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Figure 3: Correlation function resulting from simulation of 100 jets from a surface 
of a cylinder with the radius 6 fm and a "thermal" background source; 2/3 af all 
pions come from jets. Plotted are cuts along g-axes. 



In Figure El we plot the correlation function resulting from a simulation with 
100 jets where only 1/3 of the pions are produced by the background. The 
dimensions of the background source (averaged over p t ) are comparable with 
the transverse size of the source of jets, while they are clearly larger than it in 
the longitudinal direction. Therefore, for the two transverse g-components the 
background just weakens the characteristic shape from jets. On the other hand, 
in qi a narrow peak at qi = results from the large size of the background. 

We checked that if the percentage of jet pions drops as low as 1/6, the 
Pt-integrated correlation function is completely determined by the background 
source and no signal of jets is seen. Note that the proportion of jet particles can 
be even lower at ALICE: at midrapidity we assume 3000 for dN/dy of charged 
particles, out of which about 90% are pions, so we obtain roughly 1350 single- 
charge pions. Since a jet produces on average 2.4 positive pions in our acceptance 
region, 100 observable jets with pseudorapidity — 1 < r\ < 1 will leave us with 120 
single-charge pions per pseudorapidity unit. That makes the ratio of pions from 
jets to the total pion number about 120/1350 ~ 1/11. (We explore this situation 
in the next paragraph.) 
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Figure 4: Correlation function resulting from simulation of 100 jets from a surface 
of a cylinder with the radius 6 fm and a "thermal" background source; 1/12 of 
all particles come from jets and a cut is imposed: K t > 1.2 GeV/c. Plotted are 
cuts along g-axes together with Gaussian fits. 



High p t - In order to eliminate the effect of the background we can focus on par- 
ticles with high transverse momenta where we expect the portion of jet-produced 
pions to be increased. In Figure |U we plot correlation function calculated for jets 
with a background source as previously, but we use only particle pairs with K t 
above 1.2 GeV in our analysis. A non-Gaussian shape, particularly in q D and 
slightly in q s , is observed, which is due to contribution from jets. 

We have also simulated jet production from a source with Gaussian transverse 
profile together with background and imposed the same contraint on K t (not 
shown). This leads to similar non-Gaussian shapes of the correlation function as 
shown in Fig. HJ A more detailed study of the shape of the correlation function 
is neccessary in order to understand differences between the two models for jet 
distribution. 



5 Conclusions 

• Correlation functions are affected if there are pions produced from fragmen- 
tation of jets. 



8 



• Correlation functions assume particularly characteristic shape if there is 
strong suppression of jets due to energy loss of the leading parton in a 
deconfined medium. 

In the latter case, jets are effectively emitted only from the surface of the fireball. 
Presumably, these characteristic shapes are dissolved in real data because there 
are much more particles stemming from thermal fireball. 

Correlation signal of individual jets is not visible if the total number of parti- 
cles in the bin is much larger than the number of particles produced by a jet. We 
can estimate the contribution to each bin of correlation function as N 2 ot (N tot is 
total number of particles). We write N tot = Nj + Nb, where Nj is the number of 
particles produced by jets and Nj, the number of particles coming from the thermal 
source. Then the number of pairs is written as (Nj + Nb) 2 = Nj 2 + 2NjNb + Nb 2 . 
In our study, Nj corresponds to signal, N 2 to background, and NjNf, is the 
correlation between particles of signal and background. If Nj <C Nb, the signal 
becomes invisible. 

• The influence of background can be limited by using only pairs with high 
transverse pair momentum K t in the analysis. 

Although we couldn't describe the correlation function obtained after such a cut- 
off by our analytical expressions, we could observe a clear deviation from Gaussian 
shape. Note that Gaussian shape would result from our background source in the 
absence of jets. A non-Gaussian shape can influence results of fitting real data 
by Gaussian parametrisation and may cause problems when interpreting such a 
fit. We have shown that jets lead to structures in the correlation functions at 
large momentum differences: beyond the main peak at q = 0. Thus it is also 
important to look at the large-q region in order to learn about the characteristics 
of the source. 
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